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Abstract 

A new version of classical S-procedure in system theory is 
proposed based on duality in the space of positive definite 
matrices and introduction of matrix Lagrange multipliers. A 
new proof and extension of the recent results [1] concern- 
ing equivalence between frequency domain inequality on fi- 
nite frequency range and constrained dissipativity property 
for linear systems is given. The results of this paper extend 
S-procedure to allow for analysis and design of robust sys- 
tems with matrix inequalities constraints. 

Keywords. Kalman-Yakubovich-Popov (KYP) lemma, 
S-procedure, frequency domain inequality, linear matrix in- 
equality, dissipativity. 

1 Introduction 

Recently a number of new tools for systems analysis and de- 
sign related to frequency domain inequalities (FDI) over a 
finite frequency range (so called Generalized KYP-lemma) 
have been developed [2-A]. It follows from the results of 
[2—4] that fulfillness of a standard FDI in a finite frequency 
range is equivalent to validity of some nonclassical linear ma- 
trix inequalities (LMI) for a pair of matrices P, Q replacing 
inequalities for a single matrix P appearing in the classi- 
cal KYP-lemma. It was shown in [1] that FDI, in turn, are 
equivalent to some time-domain inequality (TDI, dissipation 
inequality [5]), valid only over a part of the system trajec- 
tories, determined by an additional integral matrix inequal- 
ity (restricted or constrained dissipativity [1]). Thus, a com- 
plete extension of the classical KYP-results on equivalence 
between FDI, TDI and LMI to the "finite-frequency" case was 
obtained. Note that the proof of equivalence between FDI and 
TDI in [1] goes along the lines of the necessity proof for the 
frequency-domain absolute stability criterion [6,7]. 

In this paper a new proof of the result of [1] is provided 
based on the losslessness result for a new version of the clas- 
sical S-procedure [8]. A new version of the S-procedure, also 
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included in the paper, deals with constraints in LMI form, or 
more generally, conic inequalities in linear spaces. 

In the next section a new S-procedure results are pre- 
sented. In Section 3 they are applied to the proof of equiv- 
alence between TDI and LMI. 

We use the following notation. The set of square inte- 
grable functions on [0, oo) is denoted by £2 [0, 00), Af^ where 
M is a matrix, stands for its transposition and complex con- 
jugate of all elements. For a square matrix M, its Hermitian 
part is defined by He(A/) := (Af + M'')/2. The interior of a 
set r2 is denoted by Int ft. 

2 Conic S-procedure 

Let X, Yi, . . . , Yyn be linear topological spaces, Gj : X ^ 



. . , TO be continuous mappings. 



, TO a convex cone Kj C Yj be 



Let for any j = 1, 
given defining inequality Gj (x) > for x G X as inclusion 
Gj{x) G Kj. Let Y* denote an dual space to Yj, i. e. a linear 
space of linear continuous functionals y* on Yj and K* C Y* 
denote an dual cone to A'^, i. e. K* = {y* G 1^* : (y*, j/j) > 
Vyj e Kj^^ where (y*, j/j) is the value of the functional y* 
at the element yj . 

Obviously, if F = x Fi x • • ■ x y„„ then Y* = 

X Yj^* X • ■ • X Y*^ is the set of all corteges (^y^^yl, . . . ,y^), 
where t/p G M^, y* is a linear functional from Y*. 

Consider the following two relations for the mappings Fq, 

Gl, . . . , Gm- 



(A) F{x) >OfoTx€ X, Gj{x) G Kj, j = l,...,m; 

rn 

(B) 3to > 0, Tj e K*: toF{x)- J2 {rj,Gjix)) > Vx G 



X. 



Obviously, validity of (B) with tq > implies (A). Indeed, if 
X G X satisfies inequalities Gj{x) G Kj, j — 1, . . . ,m, then 
it follows from (B) that tqF{x) > 0, since (r^, Gi{x) > O) 
for j = 1, . . . , m. The opposite statement is not true even 



in the case of scalar constraints K, = 



J = 1,- 



corresponding to the classical S'-procedure [8]. 

Similarly to the classical case we will say that S'-procedure 
with conic constraints Gj{x) > is lossless, if (B) with 
To > implies (A). 
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It is well known [9] that losslessness of the classical S- 
procedure is equivalent to the duality theorem in the corre- 
sponding optimization problem. However, the problem is, in 
general, nonconvex and only a few classes of functionals F, 
Gi, . . . , Gm are known to possess the losslessness property. 

For example, classical S'-procedure is lossless, if m = 1 
and F, Gi are quadratic forms on real or complex linear space 
X. It is also lossless, \f m — 2 and F, Gi, G2 are quadratic 
(Hermitian) forms on the complex linear space X. However, 
classical 5-procedure for quadratic forms is, in general, lossy 
for TO > 2 in real case and for m > 3 in complex case [9]. 
A. Megretski and S. Treil proved in 1990 [11] that the classi- 
cal S'-procedure is lossless for all m > 1, if F, Gi, . . . , Gm 
are integral quadratic forms on £2 (0, 00). V. Yakubovich ex- 
tended this result to a more broad class of quadratic function- 
als, forming the so-called S-system [10]. 

Below an extension of the results of [10] to the case of 
the S'-procedure with conic constraints is formulated. Note 
that the general formulation of the S-procedure with conic 
constraints was presented, e.g. in [12, 13], 

Theorem 1 . Let — 1 and K is the closure of the cone 
K generated by the set 

T{X)^{{Fo{x),Fi{x),...,Fm{x)) : x e X}. 

If the cone K is convex, then the S-procedure with conic 
constraints is lossless. 

If in addition, constraints Gj (x) G Kj are regular, namely 
3xo : Gj (xo) G Int Kj, then one can choose tq ~ I in (B). 

Proof. Condition (A) implies that F{x) > for Gj{x) € 
IntKj, i.e. intersection of the set T{X) and the open cone 
D = {(-2/0- yi, ■ ■ ■ ,ym) ■■ yo > 0,yj e IntKj.J 1, . . . ,to} 
is empty: D f] T{X) = cj). Therefore, D(^K ^ cf). Apply- 
ing separation theorem for cones, we obtain that there exists 
vector T* = (to , Tj* , . . . , r*„) e Y* such that < t^,F{x) > 
+ < T-],G{x) >> for all a; e X and < r*, y >< 

for all y^D, i.e. < r*, yo > + E7=i < ^hVj ><^- 
any j ~ 1 .... to pick up yj G liAKj and choose sequences 
2/Ofc ^ 0, ysfc -> as /c ^ 00, such that y^^ > 0,ys^ G 
Kg, s ^ j. If k 00, then we obtain < t* , yj >< 0, i.e. 
— r* G K* . The first part of the theorem is proved. 

Taking x — xq from regularity condition and yo ^ 
yields Tq , ?/o > 0. i-e. Tq > 0. Dividing the inequality (B) by 
To, we arrive at the second statement of the theorem. End of 
the proof. 

Our next step is to extend the definition of S-system [10] 
to the case of conic constraints. 

Definition 1 . Let Fj, j — 0, 1, . . . , m be mappings from 
a Hilbert space Z to spaces of self-adjoint operators over 
corresponding Euclidean space IR"^, such that Fj : 
SR{nj X rij). 

We say that Fq, Fi, .. ., F„ifonn a S-system if there exists 
a subspace Zq and a sequence of linear bounded operators 
Tk : Z Z, k — 1, 2, . . . such that 



(i) < TkZi, Z2 >— > as k 00 for all z\, Z2 G Z; 
{ ii) Zq is invariant for Tk for all k — 1, 2, . . . ; 
(Hi) Fj(Tkz) —f Fj{z) as k 00 for all j — 0, I, ...m, 

Z G Zq. 

Lemma 1 Let Fj, j — 0, 1, mform S-system. Define 

n 

the map J- : Z ^ Yi -H^"^ ^"^ nieans of the relation 
j=o 

T{z) = Fo(z) ®Fx{z)®---® Fmiz) 

Then the closure of the image .F(Z) is a convex set in 

In the special case Uj — 1, j ~ 1, . . . , to. Lemma Q co- 
incides with Lemma 1 of the paper [10] and it is proved sim- 
ilarly to the Lemma 1 of [10]. 

Example 1. An important series of examples for S-systems 
is provided by finite family of integral quadratic operators 
on the Hilbert space £2(0, cxd) of square integrable functions 
with values z{t) G K"j . The mappings are defined for any 
z G C2[0, 00) as follows: 

00 

Fjiz)=Ue J Fjz{t)z^it)Fpdt, (1) 


where F^ , F^ are Uj x symmetric matrices. In this case 
the family of the operators Tk can be chosen as time shifts: 
Tk{z){t) — z{t + k), while the subspace Zq can be chosen as 
the set of functions with zero initial conditions: 

Zo - {z(-) : z(.) G /:2(0,oo),Z2(0) = O}. 

The proof of the S-system property for Example 1 is again 
similar to [10]. Note that the cone of positive semidefinite 
matrices is selfdual. Therefore S-procedure with conic con- 
straints determined by functions ([0 deals with positive semidef- 
inite matrix Lagrange multipliers. 

Properties of the S-procedure in general case are given by 
the following theorem. 

Theorem 2 . S-procedure with conic constraints is lossless 
for any family of self-adjoint operators Fq, Fi, ... , Fm form- 
ing an S-system. 

Proof follows immediately from Theorem [J and Lemma 
[2 The result can be extended to the case of equality con- 
straints and to the case of the so called generalized S-procedure 
introduced in [2]. 
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3 Constrained dissipativity 

In this section, we first present a special case of the general- 
ized KYP lemma [4], characterizing FDIs in the continuous- 
time setting. Let complex matrices A, B, 11, and real scalars 
TUi, ZU2 be given. Define 



fl -.^ { uj e JR \ (iiJ ~ n7i)(w - VJ2) < }. 
(We may assume uj2 > without loss of generality). 



(2) 



Theorem 3 [4]. Suppose n is Hermitian matrix, pair (A, B) 
is controllable, and has a nonempty interior. Then the fol- 
lowing statements are equivalent. 



(i) The frequency domain inequality 



{jLuI - A)- 

I 



n 



{jLuI - A)- 

I 



'B 



< 



(3) 



holds for all Lu £ such that det{jLuI — A) ^ 0. 

( ii) There exist Hermitian matrices P and Q such that 
Q > and the linear matrix inequality 



'A 


B' 


* - 


I 








p 



-Q P + jn^o' 



A B 
I 



-n < 

(4) 



holds, where tUq ;= (v3\ + •n72)/2. 



holds for all solutions o/0 with u G /32[0, 00) such 
that x{0) — 0, X £ £2(0, cxd) and 

00 

He J {wix + jx){vj2X + jx)*dt < 0. (7) 



Note that the corresponding result of [ 1 ] was obtained un- 
der additional condition of asymptotic stability for the system 
(|5} which is not required in the current statement. 

In Theorem @] a general frequency interval O is consid- 
ered for the FDI, and this has translated to the input constraint 
described by Q. Though the physical meaning of this con- 
straint may be not clear in general, it becomes clear for the 
following special case. 

Corollary 1 Let real matrices A, B, 11, and a positive scalar 
ZD be given. Suppose 11 is symmetric and consider the system 
where x{t) £ IR" is the state and u{t) £ IR™ is the in- 
put. Assume that {A, B) is controllable. Then the following 
statements are equivalent. 

( i) The frequency domain inequality holds for all uj 
such that \lo\ < vd. 

( ii) The time domain inequality ^ holds for allu £ [0, cxa) 
such that X £ £2(0, 00) and 



Choosing the parameters wi = W2 — and t — the 
set n becomes the entire real numbers, and thus statement (i) 
becomes the FDI for all frequencies. In this case, the term 
associated with Q in the LMI Q becomes positive semidef- 
inite, and hence the best choice of Q for satisfaction of (@J 
is Q = 0- The resulting LMI with variable P is exactly the 
same as the one in the standard KYP lemma. 

The following result extends the result of [1]. It provides 
an equivalence between FDI and time domain dissipation in- 
equality over a restricted class of input signals. 

Theorem 4 . Let complex matrices A, B, II, and real scalars 
rui, W2 be given and fl be defined by (|2j. Consider the sys- 
tem 

x{t) ^ Ax{t) + Bu{t), tG[0,oo), (5) 

where x{t) £ (D" is the state and u{t) £ (D'" is the input. 
Assume that (A, B) is controllable, II is Hermitian, and Jl 
has a nonempty interior. Then the following statements are 
equivalent. 

(i) The frequency domain inequality (0 holds for uj £ CI. 

( ii) The time domain inequality 



f 


X 


n 


X 




u 




u 



dt<0 



(6) 



xi^dt<..^Jxx^dt. 




(8) 



Moreover, the above two statements are equivalent when the 
two inequalities "<" are replaced by ">." 

Loosely speaking, the first part of Corollary states that 
the FDI in the low finite frequency range means that the sys- 
tem possesses the property (|6j for the input signals u that 
drive the states not too fast (slowly). The bound on the "slow- 
ness" is given by w in the sense of (|8}. The second part of 
Corollary[2makes a similar statement for the FDI in the high 
frequency range. To derive Corollary ^ from Theoreml^one 
needs to put uJi ~ —u)2 — vo. 

Proof of Theorem In view of Theorem|3]it is sufficient 
to prove equivalence of (ii) to the condition (ii) of Theorem 0] 
(solvability of matrix inequality (|4}). The result follows from 
Theorem 121 of the previous section with m = 2 (one matrix 
constraint). Denote 



dt, 







00 




t 




z — 


X 


-^'"-/ 


X 


n 


X 




u 


u 




u 



00 

Gi{z) — —Yie j {vjix ^ jx){vj2X -\- jx)^ dt. 



(9) 
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Obviously, TDI (|6|l,Q correspond to the statement (A) of [5] 
the conic S-procedure. At the same time the statement (B) 
means existence of a n x n-matrix r* from the dual cone 
K* to the cone of positive semidefinite matrices satisfying 
inequaUty < r*, G'i(z) >> Vz or 



[6] 



He 



X 

u 



it 



n 



X 

u 



— {tuix + ]x)* T{m2X + ]x) \ dt >Q. 



(10) 

Replacement of r by Q and substitution of x from (|5j trans- 
forms ( I10> into LMI 0. To verify regularity condition of 
Theorem|2]take xq — {x, u), where 

x{t) = -{A + ^iy^Bexp{~fit), u{t) = exp(-/it), ^ > 

and /I G IntJl. Application of Theorem|2lends the proof. 
Remark. Similar results hold for discrete-time case. 



4 Conclusions 

The property of the system defined by the item (ii) of Theo- 
rem 0]and the Corollary can be called constrained dissipativ- 
ity or restricted dissipativity. It is weaker than standard pas- 
sivity or dissipativity conditions and may better reflect spec- 
ifications for real systems. At the same time the property 
of "slowness" described by the inequality (|8} leaves enough 
flexibility to be useful for robustness analysis of systems. 

The results of the paper shed new light on the intimate 
interrelations between S-procedure and KYP-lemma. They 
allow to extend classical S-procedure tool to allow for analy- 
sis and design of robust systems with matrix inequalities con- 
straints. 
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